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Abstract. The authors investigate the spontaneous emission dynamics of self- 
assembled InGaAs quantum dots embedded in GaAs photonic crystal waveguides. For 
an ensemble of dots coupled to guided modes in the waveguide we report spatially, 
spectrally, and time-resolved photoluminescence measurements, detecting normal to 
the plane of the photonic crystal. For quantum dots emitting in resonance with the 
waveguide mode, a ~ 21 x enhancement of photoluminescence intensity is observed as 
compared to dots in the unprocessed region of the wafer. This enhancement can be 
traced back to the Purcell enhanced emission of quantum dots into leaky and guided 
modes of the waveguide with moderate Purcell factors up to <~ 4x. Emission into 
guided modes is shown to be efficiently scattered out of the waveguide within a few 
microns, contributing to the out-of-plane emission and allowing the use of photonic 
crystal waveguides as broadband, efficiency-enhancing structures for surface-emitting 
diodes or single photon sources. 
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Photonic crystal waveguides (PCWs) are of strong interest as optical elements for 
integrated nanophotonic optical circuits and on-chip quantum optics applications. [1, 
2, 3, 4, 5] They have been used to route single photons from cavity-coupled [2, 6, 7, 5] 
and waveguide-coupled quantum dots (QDs) [8, 9], but also to tailor the local density 
of optical states (LDOS) an emitter experiences. This method to modify the LDOS 
experienced by an emitter provides a route to engineer the rate and directionality of 
spontaneous emission. [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 5, 22] This is a 
key concept to enhance the efficiency of nanoscale light sources such as single photon 
sources [23, 24, 25, 26, 27, 28] and nanoscale lasers. [29] Recently, Kaniber et al. [30] 
demonstrated a ~ 16 x enhanced extraction efficiency for single QDs emitting into 
the two-dimensional photonic bandgap of a photonic crystal. The photonic bandgap 
inhibits photon emission into the in-plane direction and redistributes it into out-of- 
plane modes, effectively increasing the extraction efficiency. However, this comes at 
the cost of long radiative lifetimes, imposing an inherent jitter in the photon emission 
time, a source of quantum distinguishability. [31, 32] In contrast, strong enhancements 
of spontaneous emission rates have been observed for QDs in low mode volume, high- 
Q cavities. [33, 34, 7] However, these systems require a sophisticated electro- [35] or 
thermo-optical [36] tuning method to spectrally bring the emitter and cavity mode into 
resonance. 

In this paper we demonstrate the advantage of using a PCW mode to enhance the 
emission rate and the extraction efficiency of an ensemble of QDs over a wide spectral 
range of ~ 18 meV in the out-of-plane direction. We compare the emission properties 
normal to the sample surface for QDs in the bulk GaAs, the photonic crystal membrane 
(PC) and the PCW region. Measurements were made for QDs in and out of spectral 
resonance with the guided PCW mode. These measurements show a strong enhancement 
of the photoluminescence (PL) intensity up to ~ 21 x for QDs spatially and spectrally 
coupled to the waveguide mode compared to QDs in the unprocessed bulk material. 
We attribute this enhancement to a combination of angular redistribution of emission 
into leaky PCW modes and Purcell enhanced emission into the guided PCW mode with 
subsequent scattering into leaky modes. 

The sample investigated was grown by molecular beam epitaxy and consists of a 
500 nm thick Alo.8Gao.2As sacrificial layer, and a 180 nm thick undoped GaAs layer 
with a single layer of nominally Ino.5Gao.5As QDs at its midpoint. The sample has a dot 
density of pqd > 50 pvaT 2 . A two-dimensional PC formed by defining a triangular array 
of air holes was realized using a combination of electron-beam lithography and reactive 
ion etching. PCWs were established by introducing line defects consisting of a single 
missing row of holes (Wl waveguide). Free standing GaAs membranes were created in 
a final wet etching step using hydrofluoric acid. 

For optical characterization the sample was mounted in a liquid He-flow cryostat 
and cooled to 10 — 15 K. For excitation we used a pulsed Ti-Sapphire laser (80 MHz 
repetition frequency, 6 ps pulse duration) tuned to the low energy absorption edge of the 
bulk GaAs (\i aS er = 815 nm). Excitation of the QDs and detection of the emitted PL 
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Figure 1. (a) and (b) Scanning electron microscope images of the investigated 
Wl waveguide - quantum dot system from the top. (c) Spatially-resolved scan 
of the photoluminescence signal performed with excitation and detection from the 
top and integrated over the 1298 - 1340 meV spectral range. The area of high 
photoluminescence intensity corresponds to the waveguide region. 

signal was done perpendicular to the sample surface using a 100 x microscope objective 
(NA=0.50). The full- width-half-maximum (FWHM) of the excitation laser spot on 
the sample was determined to be FWHM ~ 1.3 /mi, while the size of the detection 
spot had a diameter of FWHM ~ 6.0 /im. The QD PL was spectrally analyzed using 
a 0.5 m imaging monochromator and detected using a Si-based, liquid nitrogen cooled 
CCD detector. For time-resolved spectroscopy we used a Si-based avalanche photodiode 
connected to the side-exit of our monochromator providing a temporal resolution of 
~ 350 ps. 

In Fig. 1(a) and (b) we present a scanning electron microscope image of a nominally 
identical sample to the one used for the PL measurements. The photonic structure has 
a total length of 45 /mi, a slab thickness of h = 180 nm, a PC lattice constant of 
a = 270 nm, and an air hole radius to lattice constant ratio of r/a = 0.34. The two 
rectangles directly next to the photonic crystal structure serve for orientation purposes 
during optical measurements and have no influence on the photonic properties of the 
waveguide. We obtained a spatially and spectrally resolved PL map of this particular 
structure by scanning the laser spot over the surface of the sample and recording a PL 
spectrum at every position in a confocal detection geometry. A typical result of this 
scan is presented in Fig. 1 (c) where we have integrated the detected PL over the spectral 
range of the QDs (1298 meV to 1340 meV) and plotted the resulting intensity in a false 
color representation. The contours of the photonic structure can be easily recognized in 
the PL map and we observe a clear luminescence enhancement on the photonic structure 
compared to the unpatterned substrate. 
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We now investigate the origin of this PL enhancement. To do this, we performed 
photonic bandstructure calculations using the software package RSoft. [37] We use the 
appropriate parameters for GaAs of ncaAs = 3.5 and the geometric parameters of 
h/a = 0.6667 and r/a = 0.34, corresponding to the investigated Wl PCW. The result 
of this simulation is plotted in Fig. 2(a) where we plot the normalized frequency of 
the photonic bands as a function of k-vector on the path from the F point to the K' 
point. [38, 39] We present the PCW modes as blue solid lines, the slab waveguide modes 
as light gray regions, and the lossy region above the light cone is shaded in dark gray. 
The region above the light cone corresponds to the energy-wavevector combinations for 
which photons are not confined to the slab by total internal reflection, i.e. they can leave 
the waveguide in vertical directions. We calculate the guided part of the lowest energy 
waveguide (0th order) mode WM1 to be at a normalized frequency of a/A = 0.262—0.274 
which corresponds to an energy of E = 1203 — 1259 meV. For small wavevectors k < 0.28 
this mode overlaps with the region above the light cone. Photons of these wavevector- 
energy combinations can, therefore, couple to free space modes and leave the sample 
in vertical directions. This unguided part of the waveguide mode spans the normalized 
frequency range of a/A = 0.275 — 0.308 (E = 1260 — 1413 meV). The second waveguide 
(1st order) mode WM2 extends from a/A = 0.285 - 0.296 (E = 1309 - 1359 meV) with 
an additional unguided region at a/A = 0.291 — 0.296 (E = 1338 — 1359 meV). The third 
waveguide (second order) mode WM3 is at much higher energies of a/A = 0.334 — 0.337 
(E = 1534 - 1548 meV) (data not shown). 

For comparison we plot two examples of PL spectra in Fig. 2(b). The blue 
line corresponds to a spectrum recorded directly on the PCW, while the black line 
corresponds to a spectrum recorded next to the PCW on the unprocessed GaAs bulk 
material (note the x5 magnified scale). We clearly observe a strong enhancement of the 
PL intensity in the spectral range between 1300 meV and 1325 meV. For these excitation 
conditions, we detect the maximum signal of ~ 78000 cts/s at an energy of ~ 1317 meV, 
~ 21 x stronger than the ~ 3700 cts/s obtained from the unprocessed bulk material. 
The intensity oscillations in the spectrum recorded on the PCW structure, most likely 
originate from Fabry-Perot resonances due to the finite length of the waveguide. [22] We 
notice that the energy range over which we observe PL enhancement coincides very well 
with the flat part of the dispersion relation of the second energy mode, but also with 
the energy range of the unguided parts of the first and the second waveguide mode as 
shown in Fig. 2(a). 

We continue investigating this system by performing time-resolved measurements 
on the Wl waveguide for 200 different energies in the range between 1239.8 meV and 
1377.6 meV. To do this we used the spectrometer as a spectral bandpass filter of width 
~ 0.25 meV and recorded multiple decay transients as a function of energy. The 
measured data was then fitted with an automated fitting algorithm, taking into account 
the instrument response function of the experimental setup. The algorithm first tries to 
fit a bi-exponential decay but automatically switches to a mono-exponential decay when 
one of the fitting parameters is regarded unrealistic. This is the case when (i) one of the 
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Figure 2. (a) Photonic bandstructure calculations for a Wl waveguide with 
r/a = 0.34 and h/a — 0.6667. The blue lines correspond to the photonic waveguide 
modes and the light gray region to the slab waveguide modes. The dark gray region 
indicates the region above the light cone, (inset) Scanning electron microscope 
image of the investigated photonic crystal waveguide structure, (b) Comparison 
of the photoluminescence signal measured on the ensemble of quantum dots in the 
unprocessed bulk material (black line) and on the ensemble of quantum dots in the 
waveguide region (blue line), (c) and (d) Fitting parameters of the decay transients 
measured on the photonic crystal waveguide as a function of energy. In (c) the 
extracted decay times and in (d) the corresponding amplitudes are plotted. The black 
circles correspond to the decay parameters obtained from mono-exponential fits, while 
the blue triangles (red squares) correspond to the fast (slow) component of the bi- 
cxponcntially fitted decay transients, (insets) Decay transients and fits at the energy 
of 1275 meV (indicated by the red dashed line) for the upper, red-framed inset and 
1318 meV (indicated by the blue dotted line) for the lower, blue-framed inset. 
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amplitudes becomes negative, (ii) the two lifetimes differ by less than 20%, (iii) one of 
the lifetimes is shorter than 50 ps, or (iv) one of the amplitudes is more than 25 x larger 
than the other one. In Fig. 2(c) the extracted decay times are shown and in Fig. 2(d) 
the corresponding amplitudes are plotted. The blue triangles (red squares) correspond 
to the fast (slow) component of the bi-exponentially fitted decay parameters. When a 
decay is fitted mono-exponentially the resulting parameters are plotted as black circles. 

The decay dynamics can predominantly be fitted with bi-exponential parameters 
close to the resonance with the maximum PL enhancement at ~ 1317 meV, and with 
mono-exponential parameters off resonance. We interpret this observation according to 
the fact that only the emission dynamics of QDs spatially located in the waveguide region 
and spectrally in resonance with the PCW mode can be enhanced by the Purcell effect. 
However, our experiment detects contributions of spatially coupled and uncoupled QDs 
since our excitation spot size is 4 — 5x larger than the PCW width. In resonance 
with the PCW mode the decay times typically exhibit values of Tf ast = 0.57±0.1 ns and 
Tsiow = 2.7±0.1 ns. The Purcell enhanced decay times are longer than the ones that have 
been observed for QDs coupled to a photonic crystal cavity mode. [33, 40, 34, 41] We 
attribute this to the combined influence of an averaging effect over QDs that are spatially 
well and badly coupled to the PCW mode, the slope of the second PCW mode in k-space 
which leads to a lower photonic density of states, and additional broadening of the mode 
due to fabrication imperfections. [13] We also notice a general energy dependence of the 
decay times. At lower energies we observe longer decay times than at higher energies. 
We relate this to the enhanced contribution of fast multi-excitonic states to the PL 
signal at higher energies. [42, 43] The amplitude of the fast decay component, plotted in 
Fig. 2(d), nicely resembles the profile of the enhanced PL spectrum in Fig. 2(b), while 
the amplitude of the slow decay component stays almost constant (note the logarithmic 
scale). 

We present two decay transients with fits as examples in the insets of Fig. 2 (c) and 
(d). The decay transient in Fig. 2(c) - inset was recorded off resonance at the energy of 
1275 meV (indicated by the red dashed line), and the transient in Fig. 2(d) - inset in 
resonance at the energy of 1318 meV (indicated by the blue dotted line). Both transients 
exhibit a bi-exponential decay, albeit with a larger amplitude and shorter decay time of 
the fast decay component of the resonant transient. This is in good agreement with the 
extracted parameters presented in Fig. 2 (c) and (d). 

At first sight, the observation of a larger emission amplitude in the direction normal 
to the photonic crystal membrane may seem counter-intuitive. In ideal structures 
we would expect the Purcell enhanced emission to be efficiently guided away from 
the excitation spot along the waveguide, resulting in a decrease of emission detected 
in the vertical direction. However, disorder is known to lead to a scattering of 
the guided light via Anderson localization and radiation normal to the surface of 
the waveguide. [44, 45, 20] Indeed, recent near field measurements of the frequency 
dependence of the PCW mode in similar systems have revealed direct evidence for 
localization and its increasing importance for slow light modes. [46] Scattering in the 
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Figure 3. (a) Spatial image of the reflection of the laser spot on the sample surface, 
(b) Spatial image of the photoluminescence signal of the waveguide mode normal to 
the sample surface (filtered with a lOnm bandpass filter at the energy of the mode). 
The dotted white lines indicate the outline of the photonic crystal membrane. 



in-plane direction is inhibited due to the photonic bandgap of the surrounding photonic 
crystal, but scattering into modes above the light cone that can readily escape from the 
slab is still possible. [47, 48, 49] Furthermore, the finite length of the PCW structure 
also enables the QDs to emit in vertical directions due to Fabry-Perot effects. [50] In 
our structures, probably the most important reason for emission in vertical directions is 
the presence of non-guided waveguide modes with the same energy that can be accessed 
by scattering induced by disorder. QDs can emit light directly into the leaky modes, or 
propagating photons from the guided mode can be scattered into these modes due to 
disorder in the crystal. [47, 15, 48, 49, 22] Similar experimental observations have been 
made been made by Stumpf et al. [15], and recently by Hoang et al. [22] and Huisman 
et al. [46] 

In order to support these suggestions we performed additional measurements with 
a Si-based, Peltier-cooled CCD camera for spatial imaging. In Fig. 3(a) we present an 
image taken of the laser excitation spot on the sample surface at an energy of 1521 meV. 
The red (blue) curve corresponds to the cross section along the x- (y-) direction at the 
position indicated by the red (blue) line. From these plots, we can extract the FWHM 
of the laser spot to be ~ 1.3 /mi in both the x- and the y-directions. Fig. 3(b) shows a 
similar image of the PL emitted by the PCW mode in the vertical direction. In order 
to record this image, we mounted a 10 nm-bandpass filter, centered at the energy of the 
strongest PL enhancement, in front of the camera. While the PL signal recorded along 
the y-axis, perpendicular to the waveguide direction, exhibits a similar width as the 
reflection of the laser spot FWHM = 1.4 /im, the signal along the waveguide direction 
(x-axis) is much broader with FWHM = 6.1 /im. This observation strongly supports 
the assertion, that light efficiently emitted into the nominally guided waveguide mode 
is scattered into vertical directions due to the processes mentioned above. Since our 
detection spot has a size of FWHM ~ 6.0 /mi, we collect most of this scattered light 
and observe, therefore, not only PL emitted into the vertical direction, but also PL 
emitted into the guided modes and subsequently scattered into the vertical direction. 

We check for Purcell-enhanced emission by performing a spatially-resolved scan 
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Figure 4. Spatially-resolved scan of the time-resolved photoluminescence signal. 
The different panels show the fitted amplitude A s i ow and decay time t s i ow of the 
slow decay component and the amplitude A f ast and decay time Tj ast of the fast decay 
component. ®, ® and © present examples of decay transients of the spatially-resolved 
scan recorded at the position ® the bulk material, ® the photonic crystal membrane, 
and © the photonic crystal waveguide. 

of the time-resolved PL signal emitted at the peak of the amplitude (1317 meV). In 
order to do this we scan the excitation spot over the surface of the sample in 1 /im- 
steps and record the decay transients at every position over a 61x16 matrix. The decay 
transients are then fitted taking into account the instrument response function of the 
experimental setup. The same fitting algorithm as before was used, which first tries to 
fit a bi-exponential decay but automatically switches to a mono-exponential decay when 
one of the fitting parameters becomes unrealistic. The result of the fitting routine is 
summarized in the different panels of Fig. 4 that show the amplitude A s i ow and decay 
time T s i ow of the slow decay component and Af ast and Tf ast of the fast decay component. 
When the decay transient is fitted mono-exponentially, the extracted amplitude and 
decay time are plotted in the panels for the slow decay component and the corresponding 
pixel in the panels for the fast decay component are blackened. 

We can clearly recognize the outline of the PCW structure from the data in Fig. 4. 
On the unprocessed bulk material the decay transients are globally well described by 
mono-exponential decays, leading to the large blackened area around the photonic 
crystal membrane structure (see e.g. position ® on Fig. 4). We present a typical example 
of one of these decay transients in the panel of Fig. 4, marked by ®. The amplitude and 
the decay time of the mono-exponential decay are fitted to be Ab u ik = 111 ± 14 Gcts/s 2 
and Tb u ik = 0.6±0.1 ns, which results in an intensity of I^uik = A bu [ k -T bu [ k = 67±14 cts/s. 
The extracted parameters of the decays at the different positions are summarized in 
Tab. 1. 

Upon moving away from the bulk material onto the photonic crystal membrane, 
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Position 


r (ns) 


A (Gcts/s 2 ) 


/ (cts/s) 


hotai (cts/s) 


® Bulk 


0.60 ±0.1 


111 ± 14 


67 ± 14 


67 ± 14 


® Photonic Crystal Membrane 


0.81 ±0.1 
3.92 ±0.1 


106 ± 12 
40 ±4 


86 ± 14 
196 ± 18 


282 ± 32 


© Photonic Crystal Waveguide 


0.42 ±0.1 
1.74 ±0.16 


2048 ± 87 
147 ± 7 


860 ± 208 
256 ± 26 


1116 ±234 



Table 1. Overview on the extracted decay times and amplitudes for ® the bulk 
material, ® the photonic crystal membrane, and © the photonic crystal waveguide. 



but not on the Wl waveguide, the decay transients become clearly bi-exponential. A 
representative transient is presented in Fig. 4 and marked by ®. Here, we observe a 
very slow decay component due to the effect of the photonic bandgap and an additional 
faster decay component probably due to QDs located in proximity to the etched holes 
of the photonic crystal. The extracted parameters are summarized in Tab. 1. We 
observe an enhancement in the total intensity of ~ 4x compared to the intensity of 
the QD ensemble in the bulk material, which we relate to a more efficient collection of 
the emitted photons due to angular redistribution of emission. [51, 30] If we take into 
account the air fill factor of 42 % of the photonic crystal structure, i.e. 42 % less QDs, 
the emission efficiency is even a factor ~ 7x higher than that observed from the bulk 
material. 

Next, we focus on the waveguide region. Here, the oscillations in Af ast in Fig. 4 
are an artifact of the slightly rotated waveguide structure and the 1 //m steps during 
the scan. We observe a highly enhanced amplitude for both the fast and slow decay 
components, as can be seen in the example transient in Fig. 4 marked by ©. The ex- 
tracted parameters are summarized in Tab. 1. With a total intensity of 1116 ±234 cts/s 
the PL signal detected from the PCW is a factor of ~ 4x stronger than on the PC 
membrane. This enhancement corresponds to the spatial influence of the waveguide 
mode. The higher photonic density of states in spatial and spectral resonance with 
the PCW mode reduces the radiative lifetime of the electron-hole pairs via the Purcell 
effect. [23] Since radiative and non-radiative recombination are competing processes in 
self-assembled QDs and typical non-radiative recombination rates are comparable to the 
radiative emission rates of QDs emitting into the photonic bandgap, [52] a reduction 
in the radiative lifetime will effectively increase the radiative quantum efficiency. As 
discussed above, the average lifetime of the QDs emitting into the photonic bandgap is 
lengthened to ~ 4 — 5 ns in our sample. This value is comparatively small when com- 
pared to literature where values of ~ 10 — 12 ns [34, 30] and up to ~ 20 ns [17, 21] have 
been reported. We attribute this difference to a distinct non-radiative recombination 
rate in this specific sample, possibly due to non-ideal growth or etching conditions. To- 
gether with a smaller effective air fill factor at the position of the PCW, the suppression 
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of non-radiative recombination due to Purcell enhanced radiative recombination can ex- 
plain the increase in intensity detected from the waveguide region. Finally, compared to 
the emission from QDs in the bulk material, the measured intensity of 1116±234 cts/s 
from the waveguide is a factor of ~ 17x higher in good agreement to the emission en- 
hancement of ~ 21 x extracted from Figs. 1 and 2. Furthermore, the measured lifetimes 
of 0.15 — 0.5 ns are slightly reduced compared to the lifetime of QDs in the bulk material 
of 0.5 — 0.6 ns, which leads to a moderate Purcell factor of 1 — 4. 

In summary, we have demonstrated a ~ 21 x higher out-of-plane emission inten- 
sity for quantum dots spectrally and spatially in resonance with the photonic crystal 
waveguide mode, as compared to quantum dots located in the unpatterned bulk mate- 
rial. By comparing the emission dynamics for the different photonic environments we 
could identify this emission enhancement as a combination of angular redistribution of 
emission and Purcell enhanced emission of quantum dot transitions into the waveguide 
mode and subsequent scattering into vertical directions. This broadband enhancement 
of the photoluminescence rate and intensity can be used for the construction of efficient 
light emitting diodes or single photon sources for optical computation applications. 
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